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ABSTRACT 

We present the results of a deep wide-field near-infrared survey of the entire Pleiades clus- 
ter recently released as part of the UKIRT Infrared Deep Sky (UKIDSS) Galactic Clusters 
Survey (GCS) Data Release 9 (DR9). We have identified a sample of ^1000 Pleiades cluster 
member candidates combining photometry in five near-infrared passbands and proper mo- 
tions derived from the multiple epochs provided by the UKIDSS GCS DR9. We also provide 
revised membership for all previously published Pleiades low-mass stars and brown dwarfs 
in the past decade recovered in the UKIDSS GCS DR9 Pleiades survey based on the new 
photometry and astrometry provided by the GCS. We find no evidence of i\~-band variability 
in the Pleiades members larger than ~ 0.08 mag. In addition, we infer a substellar binary 
frequency of 22-31% in the 0.075-0.03 M Q range for separations less than ^100 au. We 
employed two independent but complementary methods to derive the cluster luminosity and 
mass functions: a probabilistic analysis and a more standard approach consisting of stricter 
astrometric and photometric cuts. We found that the resulting luminosity and mass functions 
obtained from both methods are very similar. We derive the Pleiades mass function in the 
0.6-0.03 M Q mass range and found that it is best reproduced by a log-normal representation 
with a mean characteristic mass of m c = 0.24+q'q3 M , in agreement with earlier studies and 
the extrapolation of the field mass function. 

Key words: Techniques: photometric — stars: low-mass, brown dwarfs; stars: luminosity 
function, mass function — galaxy: open clusters and associations: individual (Pleiades) — 
infrared: stars 



1 INTRODUCTION 

Over the past two decades star forming regions and rich, young 
open clusters have been the focal points of numerous searches for 
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2007). Among the key milestones in the study of substellar ob- 



is the study of the Initial Mass Function (IM F; Salpetei 195 
lMiller&Scaiolll979l: IScaldll986l : iKroup j |2002| ; IChabrienl2005l 
and its possible universality through the imaging and spectroscopic 
confirmation of candidate members in a wide variety of environ- 
ments. 

The UKIRT I nfrared Deep Sky Survey (UKIDSS; 
lLawrence et all [20070 is a deep large-scale infrared survey. 
All photometric observations, obtained with the Wide field 



jects, we should point out the discovery of a first brown dwarf in 
a star-forming region (o Qph; ILuhman et alJI 19971) . the identifica- 
tion of young L dwarfs dMartm et al.lll998|). and the first young T 
dwarf member of a Orionis (SOri 70; Zap atero Osorio et al. 2002). 
The common scientific driver of these surveys in young regions 



CAMera (WFCAM; Casali et al. 2007) equipped with five infrared 
filters (ZYJHK: fHewett et alj|200d) . are pipeline-processed at 
the Cambridge Astronomical Survey Unit (CASU; Irwin et al. 
2007, Irwin et al. in preparation^ The processed data are then 
archived in Edinburgh and released to the user community through 



* Based on observations made with the United Kingdom Infrared Tele- 
scope, operated by the Joint Astronomy Centre on behalf of the U.K. Parti- 
cle Physics and Astronomy Research Council, 
f E-mail: nlodieu@iac.es 



1 The survey is described at www.ukidss.org 

2 The CASU WFCAM-dedicated webpage can be found at 
http://apml5.ast.cam.ac.uk/wfcam 
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the WFCAM Science Archive (WSA; lHamblv et al.l l2008pl One 
of its components, the Galactic Clusters Survey (hereafter GCS) 
aims at covering ~1000 square degrees in 10 star-forming regions 
and open clusters down to 0.03-0.01 Mq (depending on the age 
and distance) to investigate the universality of the initial mass 
function. Each cluster is covered in ZYJHK with a second 
epoch in K to provide proper motions with accuracies of a few 
milli-arcsec per year (mas/yr). The latest GCS data release (DR) 
to date, DR9 on 25 October 2011, provides full coverage of the 
Pleiades in those five filters along with proper motions. This paper 
is the first of a series dedicated to the astrometric and photometric 
mass functions in young and intermediate-age open clusters as 
well as star-forming regions to address the fundamental issue of 
the universality of the IMF in an homogeneous manner. 

The rich Pleiades cluster has been subjected to a particu- 
larly high degree of scrutiny for several reasons. First, its mem- 
bers share a significant common proper motion compared to neigh- 
bouring stars estimated to (n a cos 5, fj-s) = (+19. 1 5, —4 5.72) 
and (+20.10 -45 .39) mas/yr by iRobichon etal] d 19991) and 
Ivan Leeuwenl d2009l) respectively, making astrometric selection rel- 
atively straightforward. Furthermore, the Pleiades is among the 
nearest well-populated open clust ers located at 1 34 pc from the 
Sun with an uncertainty of 5 pc I Johnsonlll957l; Gatewood et al.l 
l200d : |Pinfieldetal.ll200d : ISouth wort h" al l 120051)" while the im- 
proved reduction of the Hipparcos data by Ivan Leeuwenl d2009l) 
suggests a distance of 120.2+1.9 pc. In this work we adopt this lat- 
ter estimate. Its age has been determined using v arious methods, in- 
cluding the Zero-Age-Main-Sequence turn-off (Mermilliod 1981) 
and the lithium depletion boundary methods dStauffer et alj 1998). 
The latter method yielded a most probable value of 125 ± 8 Myr. 
Moreover, reddening along the line of sight to the cluster is gen- 
erally low, E(B -V) = 0.03 dO'dell. Hendry & Collier Cameron! 
1 19941) . Finally, the number of Pleiades members is large due to the 
numerous surveys dedicated to the stellar and substellar compo- 
nents of the c l uster djameson & Sldllertl 19891; lHamblv et alJl 19931: 
Staufferet ail 1 19941 ; IZapatero Osorio et all 1 19971; IStauffer et"al 
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Dobbie et al. 2002; Moraux et al. 2003; Deacon & Hambh 
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In this paper we present the Pleiades mass function derived 
from ~80 square degrees surveyed in 5 passbands, one at 2 epochs, 
by the UKIDSS GCS to provide photometry and astrometry for 
about one million sources. These data come from the latest data 
release of the UKIDSS GCS, DR9. In Section H we present the 
photometric and astrometric dataset employed to extract Pleiades 
member candidates. In Section [3] we review the list of previously 
published members recovered by the UKIDSS GCS DR9 and re- 
vise their membership. In Section [4] we outline two methods for 
deriving the cluster luminosity function. One method relies on a 
relatively conservative photometric selection followed by the cal- 
culation of formal membership probabilities based on objects posi- 
tions in the proper motion vector point diagram (Section [4. It . The 
second method applies a rigorous astrometric selection based on 
the formal errors on the proper motions for each photometric can- 
didate compared to the mean of the cluster (Section |4.2) after more 



3 The WFCAM Science Archive is accessible at the URL 
http://surveys.roe.ac.uk/wsa 
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Figure 1. The coverage in the Pleiades as released by the UKIDSS GCS 
DR9 (grey region). The holes are due to frames removed from the GCS re- 
lease due to quality control issues. Overplotted are previously known mem- 
ber candidates recovered by the GCS DR9 (rilled black dots). 

stringent multicolour cuts. In Section [5] we discuss the photomet- 
ric binary frequency in the substellar regime and compare it with 
previous estimates in the Pleiades and for ultracool field dwarfs. 
In Section [6] we discuss the A"-band variability of Pleiades cluster 
member candidates. In Section|7]we derive the cluster luminosity 
and (system) mass function and compare it to earlier estimates for 
this cluster and others, along with that of the field population. 



2 THE ASTROMETRIC AND PHOTOMETRIC DATASET 
2.1 Photometry 

We selected point sources in the full Pleiades cluster, in the re- 
gion defined by RA=50-64 degrees and dec= 18-30 degrees. We 
retrieved the catalogue using a Structure Query Lang uage (SQL) 
query similar to our earlier study of the Pleiades dLodieu et al.l 
120071) and applied it to the DR9 release of the GCS (see also the 
SQL query in appendix A o f lLodieu et all d2007h ). 

Briefly, we selected only good quality point sources in all 
passbands and included Z and Y non detections to allow for cooler 
brown dwarf candidates to be extracted. We did not impose a de- 
tection in the second /f-band either because the proper motions are 
not only computed with the two A'-band epochs but with any filter 
observed at a different time (see Section |2~2j . We limited our se- 
lection to sources fainter than Z = 11.3, Y = 11.5, J = 11.0, H = 
11.3, and Ki = 9.9 mag where the GCS saturates. The complete- 
ness limits, taken as the magnitude where the straight line fitting 
the shape of the number of sources as a function of magnitudes falls 
off, are Z = 20.2, Y = 19.8, J = 19.3, H = 18.4, Kl = 17.9, and 
K2 = 18.1 mag (Fig. O. The SQL query inclu des a cross-match 
with the Two Micron All Sky Survey (2MASS: ICutri et al.ll2003t 
ISkrutskie et al .12006 1 if available. The main change here compared 
to our previous study of the Pleiades made with the GCS DR1, apart 
from the significant increase in areal coverage, is the inclusion of 
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Figure 2. Completeness of the GCS DR9 dataset in the Pleiades cluster in 
each of the six filters. The polynomial fit of order 2 is shown as a red and 
defines the 100% completeness limit of the GCS DR9 in each passband. 



the proper motions determined from the multiple passband cover- 
age taken at different epochs by the GCS and released in DR9 (see 
Section ETH for more details). The query returned a total of 937,723 
sources over ~80 square degrees. The full coverage is displayed 
in Fig. [T] and the resulting (Z — J,Z) colour-magnitude diagram 
is shown in Fig. [3] We note that the theoretical isochrones plotted 
in the different figures of the paper were specifically computed for 
the WFCAM set of filters and kindly provided by I. Baraffe and F. 
Allard. 



2.2 Astrometry 

Proper motion measurements are available in the WFCAM Sci- 
ence Archive for UKIDSS data releases from DR9 for all the 
wide/shallow surveys with multiple epoch coverage in each field 
(i.e. the LAS, GCS and GPS). Details of the procedure will appear 
elsewhere (Collins & Hambly 2012), but here we give a brief de- 
scription of the process. 

When complete, each field imaged in the UKIDSS surveys 
is covered by a set of detector frames in various passbands with 
one passband revisited at least once. In general, these frames may 
have been taken at any time over the lifetime of the survey (at the 
time of writing, ~ 6 years) resulting in multiple epoch coverage 
for all sources. The approach taken in the WSA for computing 
proper motions is simplqj: first a set of local plane coordinates was 
established for all detections in all frames that have any system- 
atic offsets in absolute positions minimised; then proper motions 
were measured from linear least-squares fits in these local plane 
coordinates as a function of time. In the first process, for each 
set of frames, one reference frame was taken to map each other 



4 we emphasise that this work was done within the WSA for all the 
UKIDSS wide/shallow surveys, and not carried out for the sake of this paper 
only 



Figure 3. (Z — J,Z) CMD for ~8() square degrees in the Pleiades extracted 
from the UKIDSS Galactic Cluster Survey Data Release 9. Previously pub- 
lished Pleiades member candidates are overplotted as filled dots. The mass 
scale is shown on the right hand side of the diagram and extends down 
to ~0.02 M q, according to the NextGen (solid line) and DUSTY (dashed 
line) models (Bara ffe et all 1998 ; Chabrier et al.l2000l) . 



'slave' frame onto that master using all available detections and 
a simple linear 'plate' model. Working in local plane coordinates, 
the model comprises 6 coefficients allowing for independent ze- 
ropoint shifts and scale changes in both coordinates, rotation and 
non-orthogonality between the coordinate axes (shear). The reason 
for applying these local models is that the absolute astrometry of 
each frame is done with respect to relatively bright 2MASS stan- 
dard stars, so the zeropoint of any proper motions derived by sim- 
ply taking these raw positions would be defined by these relatively 
nearby stars, and would be subject to any bulk motions and/or drift 
of such stars as seen from our vantage point in space. The idea was 
to define a zeropoint for the proper motions across all surveys that 
is as close as possible to a true zero, i.e. one in which the average 
galaxy and quasar proper motions would be zero. Clearly this is not 
possible in low-latitude sightlines and/or for relatively shallow de- 
tection lists containing few identified extragalactic sources, so the 
best that can be done is to use the faintest and hence on average 
most distant stars possible. Since the number counts are dominated 
by the fainter stars, all stars were simply used. Note that rather than 
weighting these local model fits by the formal errors on each de- 
tection, unit-weight fits were chosen. Again, this is because the 
brighter, and on average nearer stars would carry the most weight 
given their low centroiding errors, but the local mapping models 
should not be biased towards the possibly drifting reference frame 
that would be defined by such stars which exhibit significant angu- 
lar motions simply because of their proximity to the Sun. 

Once a set of mapped local plane coordinates was set for every 
source paired across the detections available from the set of frames, 
a weighted linear least-squares fit is done for each coordinate as 
a function of time, resulting in four astrometric parameters (coor- 
dinates at a reference epoch along with proper motions in those 



© 2005 RAS, MNRAS 000,[T|-?? 



4 N. Lodieu et al. 



coordinates) plus formal errors and a standard goodness-of-fit pa- 
rameter (a reduced chi-squared statistic). 



3 CROSS-MATCH WITH PREVIOUS SURVEYS 

We compiled a list of Pleiades member candidates published over 
the past decades by various groups (Table [TJ to update their mem- 
bership status with the photometry and astrometry in DR9 of the 
GCS (Table IaTT ). This list will serve as starting point to identify 
new Pleiades members in the GCS data and derive the cluster lumi- 
nosity and mass functions. 

For the brightest mem bers, we used the ex tensive compila- 
tion of 1417 sources from IStauffer et all (2007) which includes 
candidate cluster members from several large-scale proper mo- 

i 1 1 — r I — 1 1 1 1 — I 

tion studies of the Pleiades dTr umpler 1921; Hertzsprung 119471; 
Artyukhina 1969; Jones 1981; Haro et al. 1982; van Leeuwen et al. 
1986l;IStauffer et al .11 199 ll: IHamblv et alj 19931; IPinfield et al.l2000l ; 
Adams et alj200ll ; |Peacon & Hambly 2004). We added candidates 
from a large number of additional papers dedicated to the Pleiades 
low-mass stars and brown dwarfs over the past 20 years, many 
sources being common to several studies which surveyed inde- 
pendently the same region of the cluster. The references are listed 
in Table QJ along with the original numbers of sources published 
by each study (All) and the corresponding numbers of candidates 
covered by the GCS (DR9). The success rate in the recovery of 
published candidate members is usually quite high (Table [TJ be- 
cause the GCS covers now the entire Plei ades cluster, as can b e 
seen from F i g. [TJ Earlier studies such as IHamblv et al.l dl993l) . 
I Adams et alj J200 llf ^ i and the list of 91 6 high -probability mem- 
ber candidate s front iDeacon & Hamblvl ((2004) are contained in 
Stauff eret alj d2007l) fnd ocused not only on Pleiades very-low- 
mass stars but also on brighter members which are saturated on 
the UKIDSS images. Thus these higher mass members are not re- 
trieved by our SQL query due to our photometric cuts at the bright 
end of th e survey. Sim i larly, man y L and T dwarf cand idates re- 
ported by iBihain et al.1 d2006l) and ICasewell etall j2007lFI are too 
faint to be detected on the GCS images. Some of the earlier Pleiades 
candidates are not recovered mainly because (Table [Jj: 

• 341 sources brighter than our saturation limits equivalent to 
58.5% of previously published sources not recovered by our SQL 
query (mainly coming from early surveys as mentioned above) 

• 185 sources (i.e. 31.7%) missing images in J, H, or Kl 

• 10 sources located in holes of the UKIDSS GCS coverage due 
to quality control issues (1.7% of all published members) 

• 8 (or 1.37%) known Pleiades candidates located beyond our 
3 arcsec cross-match limit, possibly because that they are higher 
proper motion non-members 

To quantify the completeness limit of the GCS Pleiades 
dataset, we have listed in Column 1 of Table [JJ (numbers in brack- 
ets) the numbers of previously-published candidates within the 
magnitude range probed by the GCS (see Section l2~H . The average 
of the percentages listed in the last column of TableQJ amounts for 
92.2% , not taking into account the sample of Zapatero Osorio et"aTJ 
dl997l) due to problems with their coordinates. We also consid- 
ered the most complete and updated sample of high-mass stars and 

5 We should mention th at this catalogue was not publish but is included in 
the lStauffer et alj ^2007) compilation 

6 Four objects have wrong coordinates, rectified in the erratum of this paper 
ICasewell et alj2010h 



low-m ass Pleiades member candidates published by Stauff er et al.1 
(2007). We recovered 94% of their sources, suggesting that overall 
our GCS sample is at least 92-94% complete. 

Table [ATI is provided in electronic format only with a total 
of 3196 known Pleiades member candidates reported in the liter- 
ature. From this sample we removed the multiple detections and 
kept all different names from earlier studies in the last column for 
future searches: we are left with 1379 Pleiades candidates. We give 
the GCS DR9 coordinates of these 1379 Pleiades candidates, the 
ZYJHKIK2 photometry, the proper motions in right ascension 
and declination with their respective errors as well as the \ 2 value 
which represents the reduced chi-squared statistic of the astromet- 
ric fit for each source. This parameter is equal to the usual chi- 
squared statistic (sum of normalised residuals) divided by the de- 
grees of freedom (i.e. number of data points minus the number of 
fitted parameters). The penultimate column supplies the member- 
ship probability for 1067 of the 1379 previously reported Pleiades 
candidates (see Section |4. 1.1 1 for the method). A total of 312 out 
of 1379 Pleiades candidates have no membership probabilities be- 
cause they are not classified as Pleiades members by the probabilis- 
tic selection (see Section |4~T1 for the photometric and astrometric 
criteria). They are divided into four groups as follows: 190 sources 
detected in ZYJHK but classified as proper motion non members, 
42 objects without Z + Y photometry, 74 objects without Z only, 
and 6 sources without Y only. The last column gives the old names 
used by earlier studies (names from different authors are separated 
by an underscore "_"). 

Previously-published Pleiades member candidates not recov- 
ered in the GCS DR9 are listed in Table IB II with their coordi- 
nates and old names from earlier studies. After removing common 
sources, we are left with 544 known Pleiades member candidates 
not in our sample. The large majority of these sources are either 
too bright or too faint to in the GCS database or photometric and/or 
astrometric non members of the Pleiades. 



4 NEW SUBSTELLAR MEMBERS IN THE PLEIADES 
4.1 Probabilistic approach 

4.1.1 Method 

In this section we outline the probabilistic approach we employed 
to select low-mass stars and brown dwarf member candidates of 
the Pleiades using photometry and astrometry f rom the UKIDSS 
GCSDR9. This method i s des cribed in detail in lPeacon & H ambly 
(2004j and lLodieu et"ai] < |2007|) . The main steps are: 

(i) Define the cluster sequence using candidates published in the 
literature within the area covered by the latest release of the GCS 

(ii) Make a conservative photometric cut in the (Z — J,Z) dia- 
gram to include known members and identify new cluster member 
candidates (dashed lines on the top-left panel in Fig. [6}. 

(iii) Analyse the vector point diagram in a probabilistic manner 
to assign a membership probability for each photometric candidate 
with a proper motion measurement (Sec tion [4. 1.21 . 

(iv) Obtain an illustrative list of high probability cluster mem- 
bers by choosing a specific threshold for the membership, chosen 
as p^0.6 here 

(v) Derive the luminosity and mass function using all candidates 
with membership probabilities without any threshold for a com- 
plete count of the membership 
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Table 1. Updated membership of Pleiades member candidates published 
in the literature and recovered by the GCS DR9. Papers dedicated 
to the Pleia des over the past tw o decades are ordered by ye a r. Ref- 
erences are: iHamblyetalJ j 19931). Zapatero Osorio et al. I (19971). I Festin 



Bouvier 



e p ast tw o 
HI [\ | 993l). 
et al l |l998l) . IZap; 
etalj j2000l).lMor; 



patero Osorio et al 
fcOOd). iMoraux et all bpolll 
Moraux et all 12003). Deaco n & Hambl 



jl999l) .IHamblv et al 
Jameson et all 1 200: 
nblvl 



Pinfield ( 
Dobbie et al l <2002l) . 

Bihain et alj J200el> lLodieu et alj <2007t) . ICasewell et alj <200j) . 
Stauffer et all i2007l) . Columns 2 and 3 give the numbers of sources pub- 
lished by the reference given in Column 1 and the numbers of sources re- 
covered in GCS DR9, respectively. Numbers in brackets in Column 1 repre- 
sent all sources in a given catalogue within the magnitude range probed by 
the GCS. Columns 4 and 5 give the numbers of high-probability members 
(p>60%) and non members (NM) according to our probabilistic approach 
(first number) and method #2 (second number). The last column gives the 
percentages of sources recovered in the GCS DR9 (i.e. Column 2 divided 
by numbers in brackets in Column 1). 



Survey 


All 


DR9 


Memb 


NM 


% 


Hamblyl993 


440 


(440) 


418 


303/302 


115/116 


95.0 


Zapatero 1997a 


9 


(9) 


2 


1/1 


1/1 


22.2 


Festinl998 


45 


(44) 


37 


24/24 


13/13 


84.1 


Staufferl998 


20 


(19) 


16 


11/10 


5/6 


84.2 


Bouvierl998 


26 


(26) 


25 


12/15 


13/10 


96.2 


Zapaterol999 


46 


(44) 


38 


8/13 


30/25 


86.4 


Hamblyl999 


9 


(9) 


9 


5/6 


4/3 


100.0 


Pinfield2000 


339 


(338) 


320 


185/187 


135/133 


94.7 


Moraux2001 


25 


(25) 


25 


12/15 


13/10 


100.0 


Dobbie2002 


90 


(87) 


61 


8/9 


53/52 


70.0 


Moraux2003 


109 


(108) 


107 


74/74 


33/33 


99.1 


Deacon2004 


916 


(746) 


674 


467/450 


207/224 


90.3 


Bihain2006 


34 


(31) 


28 


11/14 


17/14 


90.3 


Lodieu2007 


456 


(456) 


454 


376/376 


78/78 


99.6 


Casewell2007 


23 


(16) 


16 


4/7 


9/6 


100.0 


Stauffer2007 


14ie 


(944) 


888 


567/639 


321/249 


94.0 



4.7.2 Membership probabilities 

In order to calculate formal membership prob abili ties we used 
the sam e technique as iDeacon & Hamblvl [2004) and lLodieu et alj 
d2007t) to fit distribution func tions to proper motion vector point 
diagrams ( IHamblv et al.lll995l) . We refer the reader to those papers 
for more details and additional equations. First we have rotated the 
vector point diagram so the cluster lies on the y-axis using the ro- 
tation transformation below (Equations Q]and[2): 



fi x i = 0.03896 Xfi x - 0.921 x fi y 
H v i = 0.03896 x n y - 0.921 x fi x 



(1) 
(2) 



corresponding a rotation angle of 23.7 degrees, assuming a rel- 
ative proper motion of (19.7,-44.82) mas/yr for the Pleiades as 
measured on the vector point diagram created from the GCS 
DR9 data (s lightly different from the Hipparcos absolute motion; 
van Leeuwen 2009). It is common to refer to the proper motions in 
the rotated vector point diagram as p x i and fi y i. 

We have assumed that there are two contributions to the total 
distribution cj>(fi x , fly), one from the cluster ((j> c (fix, fi y )) and one 
from the field stars ((f)f([i x , f-i y )). The fitting region was delineated 
by —50 < fix < 50 mas/yr and 20 < fi y < 70 mas/yr. These were 
added by means of a field star fraction / to yield an expression for 
<j) given in Equation[3] 
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Figure 4. Vector point diagram showing the proper motion in right ascen- 
sion (x-axis) and declination (y-axis) for previously-known member candi- 
dates recovered by the GCS DR9 (filled black dots) and the new member 
candidates with membership probabilities higher than 60% identified in this 
work (red dots). 



We have assumed that the cluster distribution is characterised 
by a bivariate Gaussian with a single standard deviation a and mean 
proper motion values in each axis fj, xc and fi yc (Equation[4): 



oc exp 



(f^x fixe) ~f~ (fly M?/ c ) 

2^2 



(4) 



4>{flx,fly) = f<t>!{^x,fl y ) + (1 - f)<f>c(flx,fly) 



(3) 



The field star distribution was fitted by a single Gaussian in 
the x axis (with standard deviation T, x and mean p x f) and a de- 
clining exponential in the y axis with a scale len gth r. The use of 
a dec lining exponential is a standard method (e.g. ljones & StauffeJ 
[l99lh and is justified in that the field star distribution is not simply a 
circularly-symmetric error distribution (i.e. capable of being mod- 
elled as a 2d Gaussian) - rather there is a prefered direction of real 
field star motions resulting in a characteristic velocity ellipsoidal 
signature, i.e. a non-Gaussian tail, in the vector point diagram. This 
is best modelled (away from the central error-dominated distribu- 
tion) as an exponential in the direction of the antapex (of the solar 
motion). 

The best fitting set of pa rameters were chosen us ing a max- 
imum likelihood method (see iDeacon & Ha mblv 2004). However 
in a deviation from this method we did not fit for the standard de- 
viation of the cluster proper motions (a). Instead we calculated the 
mean astrometric error for all objects in each magnitude range and 
used this as our cluster standard de viation. This fitting process was 
tested bv lDeacon & Hamblvl d2004l) where simulated data sets were 
created and run through the fitting process to recover the input pa- 
rameters. These tests produced no signifi cant offsets in the parame- 
ter values (see Table 3 and Appendix A of IDeacon & Hamblvl2004l 
for results and more details on the procedure). Hence, we have cal- 
culated the formal membership probabilities as, 



(5) 



f<t>f + (1 " f)<t>c 
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Table 2. Summary of the results after running the programme to derive 
membership probabilities. For each Z magnitude range, we list the number 
of stars used in the fit (Nb), the field star fraction f, and parameters describ- 
ing the cluster and field star distribution. Units are in mas/yr except for the 
number of stars and the field star fraction f, The cluster star distribution is 
described by the mean proper motions in the x and y directions (p Xc and 
li yc ) and a standard deviation a. Similarly, the field star distribution is char- 
acterised by a scale length for the y axis (t), a standard deviation T, x , and 
a mean proper motion in the x direction (/z^ , ). 



z 


Nb 


f 


a 






T 






12-13 


486 


0.86 


3.22 


-1.36 


44.65 


15.51 


17.97 


-9.48 


13-14 


853 


0.79 


3.14 


0.21 


46.25 


14.95 


17.46 


-9.22 


14-15 


1394 


0.78 


3.21 


-0.11 


46.12 


13.50 


16.95 


-9.13 


15-16 


1712 


0.86 


3.28 


0.47 


45.60 


13.19 


16.06 


-9.45 


16-17 


1885 


0.95 


3.41 


0.68 


45.25 


11.16 


15.16 


-9.79 


17-18 


1256 


0.96 


3.70 


-0.26 


44.99 


9.58 


14.28 


-10.34 


18-19 


493 


0.93 


4.28 


-1.23 


44.07 


11.02 


15.09 


-11.18 


19-20 


358 


0.92 


5.93 


-0.63 


45.76 


10.30 


15.72 


-9.86 


20-21 


369 


0.84 


9.62 


7.80 


41.75 


10.80 


17.49 


-7.76 



mainly because candidates identified in earlier surveys focusing on 
Pleiades brown dwarfs did not have as much as information as the 
GCS i.e. 6-band photometry and accurate astrometry down to ~25 
Jupiter masses. Hence, our new dataset allows us to reject many of 
the earlier substellar candidates either on photometric or astromet- 
ric grounds. 



4.2 Photometry and proper motion selection 

In this section we outline a more widely used method (referred to 
as method #2 in the rest of the paper) that we applied to select low- 
mass and substellar Pleiades member candidates. This procedure 
consists of selecting cluster candidates by applying stricter pho- 
tometric cuts in various colour-ma gnitude diagrams supplem ented 
by a proper motion selection (e.g. lLodieu et ahlboOfj, 12007). This 
method is complementary but independent from the probabilistic 
approach presented in the previous section. 

First, we plotted several colour-magnitude diagrams (Fig. [5} 
to study the position of known Pleiades members identified in ear- 
lier studies and published over the past decades (Table [T}. Based 
on these known members, we define a series of lines to select pho- 
tometrically member candidates with photometry in ZYJHK in 
four colour-magnitude diagrams as indicated below (dotted lines in 
the diagrams in Fig. [6): 



• 


(Z 


- J,Z) = 


(0.50,12.0) to (1.05,17.0) 


• 


(Z 


- J,Z) = 


(1.05,17.0) to (2.40,21.5) 


• 


(Z 


- K,Z) = 


= (1.20,11.5) to (1.95,17.0) 


• 


(Z 


- K,Z) = 


= (1.95,17.0) to (4.00,21.5) 


• 


(Y 


- J,Y) = 


(0.30,11.5) to (0.55,16.5) 


• 


(Y 


- J,Y) = 


(0.55,16.0) to (1.40,20.5) 


• 


(J 


- K,J) = 


(0.70,11.0) to (0.70,16.5) 


• 


(J 


- K,J) = 


(0.70,16.5) to (1.70,19.0) 



These photometric cuts remain conservative and the contami- 
nation to the blue side of the Pleiades sequence is still high. Hence, 
the second step consisted of applying a proper motion selection in 
the vector point diagram (Fig. [4) to improve on the photometric se- 
lection. We applied a 3<r selection given the formal errors on the 
individual proper motions for each object, implying a complete- 
ness higher than 99% for normally-distributed errors. We assumed 
a mean proper motion of (+18, —42) ma s/yr for the Pleiades , 
slight ly different fro m the Hipparcos values dRobichon etalj 19991: 
van Leeuweni r2009) because the proper motion measurements are 
on a relative system rather than the Hipparcos system as described 
in Collins & Hambly (2012). The main advantage of this method is 
that it doesn't rely on a single radius for the proper motion selec- 
tion but rather takes into account the increasing uncertainty on the 
proper motion measurements for fainter stars, and allows for differ- 
ent time baselines of the epoch frames affecting the proper motion 
errors. 

After applying both the photometric and proper motion se- 
lections, we found a large number of objects lying to the blue of 
the Pleiades sequence. Hence, we applied an additional photomet- 
ric cut in the (Z — J,Z) colour-magnitude diagram, eliminating all 
sources in the Z = 12-18 mag range and located to the left of a 
line defined by (Z - J,Z) = (0.6,12.0) to (1.2,16.5). This selec- 
tion yielded a total 1 147 low-mass stars and brown dwarfs with Z 
magnitude ranging from 12 to 21.5 (Table ICTV This total number 
is similar to the number of high probability Pleiades member can- 
didates identified via the probabilistic approach; we note however 
that the membership count from the different methods is identical 



As the astrometric errors are a function of magnitude we split 
the sample into nine bins. Each bin was one magnitude wide and 
the constituent stars used to fit for six of the seven par ameters in 
the same way as described in iDeacon & Hamblvl [2004). We used 
bins of one magnitude to have a sufficient number of Pleiades mem- 
bers in each range. In the faintest bins where the astrometric errors 
increase rapidly, the number of cluster stars was so small that we 
fixed the location of the cluster on the vector point diagram (/ji C 
and jj, yc ) to the values from a brighter bin. The other parameters 
were fitted as normal. A summary of the fitted parameters from the 
prob abilistic analysis described above is given in Table [2] 

Deacon & Hamblvl ^20041) test the reliability of this method in 
Section 2.6 and find that the method accurately recovers the pa- 
rameters from simulated datasets. We refer the reader to Table 3 in 
this paper for approximate errors on parameters but note that these 
tests were performed on simulated datasets with a much larger frac- 
tion of cluster stars. Hence our cluster parameters will likely have 
higher errors than those quoted. Some parameters such as r and a x 
vary between magnitude bins due to noise as a result of two com- 
peting factors. As we go fainter the real proper motions of our field 
star contaminants decrease as our typical field star will be more 
distant but the measurement errors on the proper motions increase. 
Hence these two parameters which describe the proper motions of 
the field stars initially fall with increasing magnitude but as the 
measurement errors (as traced by the parameter a) blow up, these 
in turn also increase. 



4.1.3 Probabilistic sample 

The probabilistic approach yielded a total sample of 8797 sources 
with membership probabilities assigned to each of them. This sam- 
ple contains 947 sources with membership probabilities higher 
than 60% listed in Table [CT] Relaxing the probability to 50% 
yields a sample of 1076 Pleiades member candidates. These high- 
probability members are displayed in Fig. |5]along with previously 
published Pleiades candidates. 

We note that the sequence of member candidates in our proba- 
bilistic sample is very similar to the sequence of previously-known 
members down to Z = 16 mag (Fig. |5j- However, we observe dif- 
ferences in the Z = 16-17 mag range and at fainter magnitudes, 
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Figure 5. Colour-magnitude diagrams showing the Pleiades member candidates previously reported in the literature (black dots) and the new ones extracted 
from our probabilistic analysis (red dots). Upper left: (Z — J,Z); Upper right: (Z — K,Z); Lower left : (Y — J,Y); Lower right: (J — K,J). Overplotted 
are the 120 Myr NextGen (solid line; lBaraffe et all 19981) and DUSTY (dashed line; IChabrier et al.l200(il) isochrones shifted to a distance of 120 pc. The mass 
scale is shown on the right hand side of the diagrams and spans ^0.6-0.02 Mq, according to the 120 Myr isochrone models. 



within the counting errors if we sum the membership probabilities 
of all stars, as expected. 



4.3 Search for lower mass members 

In this section we search for fainter and cool substellar members 
of the Pleiades by dropping the constraint on the Z-band detection 
and later the Z + Y bands. 



4.3.1 Y JHK detections 

To extend the Pleiades cluster sequence to fainter brown dwarfs and 
cooler temperatures, we searched for potential candidate members 
undetected in Z. We imposed similar photometric and astrometric 
criteria as those detailed in Section |4~2l but imposed a Z non detec- 
tion and associated criteria as described below: 



No Z detection 

Y 5? 18 and J < 19.3 mag 
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Figure 6. Same as Fig.[5]but only for Pleiades member candidates selected using method #2. The YJHK and JHK-only detections have been added too. Our 
photometric criteria listed in Section l4T2l are represented by dotted lines. The faint YJHK and JHK detections are highlighted with triangles and diamonds 
surrounding the black dots, respectively. 



• Candidates should lie above the line defined by (Y — J,Y) = 
(0.55,16.0) and (1.40,20.5) 

• Candidates should lie above the line defined by (J — K,J) = 
(0.75,16.5) and (1.70,19.0) 

• The position on the proper motion vector point diagram of 
each candidate should not deviate from the assumed cluster proper 
motion by more than 3a 

This selection returned 22 additional Pleiades member candidates 
which have been added to Table ICTI along with sources identified 
with both selection methods presented earlier. All but four of them 



are indeed invisible in the Z-band images and look well detected in 
the other bands after checking the GCS DR9 images by eye. 

4.3.2 JHK detections 

We repeated the procedure described above looking for Z and Y 
non detections. We applied the following criteria: 

• No Z and Y detection 

• J = 18-19.3 mag 

• Candidates should lie above the line defined by (J — K,J) = 
(0.75,16.5) and (1.70,19.0) 
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• The position on the proper motion vector point diagram of 
each candidate should not deviate from the assumed cluster proper 
motion by more than 3a 

This query returned 19 new Pleiades candidate members. Af- 
ter checking the images by eye, we kept only one of them because 
one is actually visible in Y (although a Y detection is not reported 
in the GCS DR9 catalogue), one is visible in both the Z and Y im- 
ages. The remaining 16 have no Z and Y images available so we 
can't confirm if they are indeed drop-outs (Table |CTV 



5 THE SUBSTELLAR BINARY FREQUENCY 

The multiplicity in the substellar domain at different ages provides 
one way to constrain the formation me chanisms of brown dwarfs. 
As in our earlier study of the Pleiades dLodieu et al.ll2007l) . we in- 
vestigated the binary frequency of Pleiades brown dwarfs using the 
photometry and colours from the GCS. However, our sample is now 
two times larger and, with the proper motions, is of much higher 
quality because it is drawn from the same homogeous survey. We 
consider in this section all Pleiades member candidates selected 
through method #2. 

W e applied the same method as described in lLodieu et al.l 
d2007h to select substellar binary candidates in the Pleiades. We 
briefly summarise the technique here. Figure|7]displays two colour- 
magnitude diagrams used to identify binary candidates because 
of the large colour range and the presence of sources above the 
cluster (single star) sequence. We started off our selection in the 
(Y — K\,KX) colour-magnitude diagram (left-hand side plot in 
Fig. |7]( because it shows large colour difference in the substellar 
regime. We applied the following method: for a given magnitude, 
e.g. K\ = 15.5-16.5 mag, we defined two horizontal lines inter- 
cepting the mean value of the single object sequence. From the in- 
tercept points we defined two vertical lines with a length of 0.75 
mag (dashed lines in Fig. [7). Then we divided the box formed by 
both sequences and the vertical lines into two boxes: single stars 
lie in the bottom part whereas binary candidates in the top one. 
Except for one system which appears to the blue of the single-star 
sequence in the (J — K\,KV) diagram, the location of the binary 
candidates is confirmed in this diagram and others as well, adding 
credence to their potential multiplicity. Our method is corrobo- 
rated by the presence of four known Pleiades bro wn dwarf bina- 
ries ( I PMBD25, IPMBD 29, PP115, CFHT-Pl-IZ4lBasri & Martfr] 
1 19991 : iMartm et al.l |2000L 120031 : iBouv et alj|2006h ab ove the sin- 
gle star sequence where Teide 1 dRebolo et al.ll 1995b . an isolated 
Pleiades brown dwarf is located (red symbol in Fig.[7). Note that we 
observe a dispersion of 0.5 mag in the single-star sequence which 
can explained by the tidal radius of the cluster leading to a variation 
of 10% in the distance of the mem bers, i.e. ~12 pc corresponding 
to ~0.2 mag jPinfield et alj2000l) . 

The binary fraction was then defined as the number of binaries 
divided by the total number of objects (single stars+binaries). We 
counted 51 binary candidates (Table iDTt and 137 single stars in the 
K\ ~ 14.32-16.27 mag range, corresponding to masses between 
0.075 and 0.03 M© at the age and distance of the Pleiades. Hence, 
we derive a binary frequency of 51/(137+51) = 27.1 ±5.8% in this 
mass range for projected separations smaller than ~ 100-200 au. 
This value is likely a lower limit because some binaries may hide 
in the single star sequence due to higher separations or mass ra- 
tios and we have not considered wider systems. The overall result 
is in agreement within the error bars wi th our previous es timate 
although on the lower side (36.5+8.0% lLodieu et al]|2007h . We 



also divided up this mass range into two bins covering 0.075-0.05 
M and 0.05-0.03 M , yielding binary fractions of 33.0±9.1% 
(35 binaries and 71 singles) and 19.5±7.0% (16 binaries and 66 
singles), respectively. The binary frequency over the lowest mass 
range is not reproduced by the latest hydrodynamical simulation of 
a 500 Mq stellar cluster as no bro wn dwarf b inary was found in 
the 0.07-0.03 M range (Table 2 in lBatj201l1) . Unfortunately, we 
cannot test further those theoretical predictions with estimates on 
the numbers of triple and high-order multiple, separation distribu- 
tions, and mass ratios. 

We investigated the range of validity of our binary frequency 
by plotting the expected positions of binary systems with primary 
masses of 0.075, 0.05, and 0.03 Mq (blue crosses in Fig. [7] Adding 
smaller mass brown dwarfs (going from 0.075 M© i.e. equal-mass 
binaries down to 0.02 Mq) to primaries with masses of 0.075 Mq 
places those systems to the red of single stars and then turn over 
towards higher luminosities. The same behaviour is observed for 
binaries with primary masses of 0.05 Mq whereas binaries with 
primary masses of 0.03 Mq turn redder and brighter. To estimate 
the sensitivity of our binary frequency in terms of mass ratios, we 
considered the photometric errors of the GCS for the three mass 
values. We inferred that our binary frequencies over the 0.075-0.05 
Mq and 0.05-0.03 M© are valid for mass ratios larger than 0.4-0.5 
and 0.8, respectively. Hence, the factor of two observed in the mul- 
tiplicity between those two mass bins should be interpreted with 
caution because they are not valid over the same mass ratio inter- 
val. If we consider only binaries with mass ratios larger than 0.8 
in the 0.075-0.05 Mq mass range, we derive a binary fraction of 
9/(71+9) = 1 1.3+5.0%, lower by a factor of two than the frequency 
inferred in the 0.05-0.03 M© range although consistent within the 
error bars. 

We note that we applied the same procedure to the probabilis- 
tic sample and found very similar numbers which do not change 
the main conclusions discussed in this Section. We infer a substel- 
lar binary frequency of 24.3+7.3% which can be divided up into 
29.2+12.0 and 20.0+8.7% for the 0.075-0.05 M© and 0.05-0.03 
M© mass bins, respectively. 

Our binary fraction is higher by a factor two to three than 
the frequency inferred from high-resolutio n imaging (13.3l"4 3 %; 
IMartm et al.M2000l 120031 : IBouv et al.| [2006) especially if we con- 
sider that all known Pleiades brown dwarf binaries lie in the 0.065- 
0.055 Mq mass r ange. Our mul t iplicit y is lower than the 50% es- 
timate derived bv lPinfield et al.l 12000) from a purely photometric 
estimate (i.e. no proper motion measurements in volved) and on the 
low s ide of Monte-Carlo predictions (32-45%; iMaxted & Jeffries! 
2005). However, our r esults are in agreement with the upper limit 
of 26+10% quoted bv lBasri & Reinersl d2006l) for low-mass stars 
and brown dwarfs. This total binary frequency is split into 11% 
of spectroscopic binaries with projected separation below 3 au, 
and 15% of wider binaries (3-15 au) in agreement with high- 
res olution imaging survey s of field ultracool dwarfs (for a review, 
see lBurgasser et a l. 2007) and theoretical predictions from hydro- 
dynamical simulations ( Batd l201 If) . 



6 VARIABILITY AT 120 MYR 

In this section we discuss the variability of the Pleiades low-mass 
stars and brown dwarfs using the two epochs provided by the GCS. 
We only considered the Pleiades member candidates with mem- 
bership probabilities larger than 60%, many of them being already 
published in the literature (Tables lATl and lCTV 
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Figure 7. (Y — K,K) and (J — K,K) colour-magnitude diagrams showing all Pleiades member candidates selected with method #2 as well as the YJHK 
and JHK-only detections. The multiple system candidates identified photometrically in the (Y — K,K) diagram are marked with a circle; the red square 
shows the location of the prototype Pleiades brown dwarf Teide 1. The mass scale is shown on the right hand side of the diagrams for an age of 120 Myr and a 
distance of 120.2 pc. Photometric error bars are shown as crosses on the right-hand side of the panel. The method describing the selection of binary candidates 
in the (Y — K,K) diagram is highlighted with black dashed lines. Blue crosses represent the expected sequences of binary systems for primaries with masses 
of 0.075, 0.05, and 0.03 M , respectively. 



Figure [8] shows the (A1-A2) vs Al for all Pleiades member 
candidates with probabilities higher than 60%. The brightening in 
the Al = 11-12 mag range is due to the difference in saturation 
characteristics between the first and second epoch, of the order of 
0.5 mag both in the saturation and completeness limit. This is un- 
derstandable because the exposure times have been doubled for the 
second epoch with relaxed constraints on the seeing requirement 
and weather conditions. We excluded those objects from our vari- 
ability study. Overall, the sequence indicates consistent photometry 
between the two K epochs and very few objects appear variable in 
the A-band. 

We selected variable objects by looking at the standard de- 
viation, defined as 1.48 x the median absolute deviation which is 
the median of the sorted set of absolute values of deviation from 
the central value of the K\-K2 colour. In the Kl = 11-13 mag, 
we identified four potential variable abjects with differences in the 
A'-band larger than 3a above the standard deviation. However, all 
appeared saturated in the second epoch images suggesting that the 
variability may be caused by the inaccurate photometry derived 
from saturated sources. We extracted another potential variable 
low-mass star around A~14.5 mag but this source is located at 
the edge of the detector, casting doubt on any intrinsic variability. 
The analysis is not possible beyond 16 mag due to the small number 
of Pleiades members with high probabilities. 

We conclude that the level of K-band variability at 120 Myr 
is small, with standard deviations in the 0.05-0.08 mag range, sug- 
gesting that it cannot account for the dispersion in the cluster se- 
quence. The same conclusions are drawn from the sample derived 
from method #2. 



7 THE INITIAL MASS FUNCTION 

In this section we discuss the cluster luminosity and mass functions 
derived from the samples of Pleiades member candidates extracted 
from both methods described in the previous section. We did not 
attempt to correct the mass function for binaries. 

7.1 The cluster luminosity function 

In this section, we construct two luminosity functions: i) we used 
the sample of 8797 Pleiades cluster member candidates selected 
by the probabilistic approach (Section [4.1) ; and ii), the 1147 can- 
didates identified with method #2 (Section 14. 2\ . The luminosity 
function of the former method is derived by summing membership 
probabilities of all stars fitted to distribution functions in the vec- 
tor point diagram, whereas the luminosity function of the latter is 
derived simply by summing the number of member candidates. 

Both luminosity functions i.e. the number of stars and brown 
dwarfs as a function of magnitude plotted per 0.5 mag bin is dis- 
played in Fig. [9] The brightest bin is a lower limit due to the satura- 
tion limit of the GCS survey. The last bin is very likely incomplete 
due to the constraint imposed on the Z-band detection. The num- 
bers of objects per 0.5 mag bin increase quickly to reach a peak 
around Z = 14.5-15 and drop off afterwards down to the complete- 
ness of our survey with a possible peak beyond Z = 20 mag (Ta- 
bles[3]&|4jl- Both luminosity functions look very similar and match 
each other within the error bars. Therefore, we conclude that both 
methods provide the same result and a good representation of the 
Pleiades luminosity and mass functions. 

Assuming that the observed lithium depletion boundary 
is at M ~ 0.075 M R QAz = 11.44 IStauffer et al] 1 19981; 
iBarrado v Navascues et al. and a distance of 120.2 pc, the 
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Figure 9. Luminosity (left) and mass (right) functions derived from our analysis of the UKIDSS GCS DR9 sample of Pleiades member candidates. The 
left-hand side panel compares the luminosity function obtained from the probabilistic approach (black symbols and black line) and the luminosity function 
derived from the selection outlined by method #2 (red colour). Error bars on the number of objects per magnitude bin are Gehrels errors. The right-hand side 
panel compares the Pleiades mass function derived from the probabilistic approach (filled black dots linked by a solid line) and the mass function derived from 
method #2 (blue symbols and blue line). The black and blue lines are log-normal fits to the Pleiades mass function. Error bars on the dN/dlogM are Gehrels 
errors whereas error bars on the masses come from the 3cr lower and upper limits on the age (100-150 Myr) and distance (114.5-125.9 pc). The first and 
last two points of the Pleiades luminosity and mass functions are to be treated with caution due to saturation and contamination at the bright and faint e nds, 
respectively. The mass function from our GCS DR1 analysis (green diamond and green line) and the field mass functions (red lines; Chabrier 2003, 2005) are 
also included for completeness and comparison. 



Table 3. Values for the luminosity and mass functions (both in linear and logarithmic scales) per magnitude and mass bin for the Pleiades open cluster using 
the probabilistic approach (method 1). We assume a distance of 120.2 pc and employed the NextGen and DUSTY 120 Myr theoretical isochrones for the mean 
values. The uncertainties given in brackets for the mid-mass (Column 4) come from the 3cr lower and upper limits on the age and distance, assuming errors of 
8 Myr and 1.9 pc, respectively. 



Mag range Nb.obj Mass range Mid-mass dN errH errL dN/dM errH errL dN/dlogM errH errL 



12.0-12.5 


179 


0.6200-0.5400 


0.5800 (0.5550-0.5890) 


16.18 


5.11 


3.99 


202.29 


63.94 


49.90 


2.43 


0.27 


0.28 


12.5-13.0 


306 


0.5400-0.4690 


0.5045 (0.4850-0.5125) 


51.03 


8.20 


7.13 


718.70 


115.43 


100.36 


2.92 


0.15 


0.15 


13.0-13.5 


372 


0.4690-0.3890 


0.4290 (0.4025-0.4380) 


81.18 


10.05 


9.00 


1014.71 


125.64 


112.45 


3.00 


0.12 


0.12 


13.5-14.0 


480 


0.3890-0.3030 


0.3460(0.3145-0.3595) 


91.95 


10.63 


9.58 


1069.15 


123.58 


111.35 


2.93 


0.11 


0.11 


14.0-14.5 


614 


0.3030-0.2300 


0.2665 (0.2365-0.2800) 


135.02 


12.65 


11.61 


1849.55 


173.31 


159.03 


3.05 


0.09 


0.09 


14.5-15.0 


779 


0.2300-0.1740 


0.2020 (0.1775-0.2140) 


158.64 


13.63 


12.59 


2832.91 


243.31 


224.74 


3.12 


0.08 


0.08 


15.0-15.5 


848 


0.1740-0.1350 


0.1545 (0.1365-0.1655) 


139.02 


12.82 


11.78 


3564.67 


328.78 


302.06 


3.10 


0.09 


0.09 


15.5-16.0 


863 


0.1350-0.1070 


0.1210(0.1073-0.1305) 


94.90 


10.78 


9.73 


3389.36 


385.01 


347.46 


2.97 


0.11 


0.11 


16.0-16.5 


931 


0.1070-0.0862 


0.0966(0.0858-0.1048) 


45.73 


7.82 


6.74 


2198.46 


375.84 


324.22 


2.69 


0.16 


0.16 


16.5-17.0 


953 


0.0862-0.0710 


0.0786 (0.0671-0.0855) 


38.41 


7.26 


6.18 


2527.30 


477.51 


406.43 


2.66 


0.17 


0.18 


17.0-17.5 


747 


0.0710-0.0549 


0.0630 (0.0534-0.0711) 


35.81 


7.05 


5.96 


2224.41 


437.69 


370.40 


2.51 


0.18 


0.18 


17.5-18.0 


508 


0.0549-0.0486 


0.0517(0.0463-0.0589) 


16.22 


5.12 


4.00 


2574.60 


812.61 


634.33 


2.49 


0.27 


0.28 


18.0-18.5 


268 


0.0486-0.0440 


0.0463 (0.0414-0.0506) 


16.70 


5.18 


4.06 


3630.65 


1125.53 


881.74 


2.59 


0.27 


0.28 


18.5-19.0 


224 


0.0440-0.0396 


0.0418 (0.0383-0.0461) 


13.66 


4.80 


3.66 


3105.68 


1090.16 


832.42 


2.47 


0.30 


0.31 


19.0-19.5 


164 


0.0396-0.0368 


0.0382 (0.0356-0.0420) 


11.34 


4.48 


3.33 


4051.07 


1599.11 


1189.51 


2.55 


0.33 


0.35 


19.5-20.0 


193 


0.0368-0.0339 


0.0353 (0.0330-0.0390) 


15.26 


5.00 


3.87 


5260.69 


1724.40 


1335.78 


2.63 


0.28 


0.29 


20.0-20.5 


223 


0.0339-0.0311 


0.0325 (0.0305-0.0370) 


38.49 


7.26 


6.18 


13745.71 


2594.29 


2208.46 


3.01 


0.17 


0.18 


20.5-21.0 


145 


0.0311-0.0294 


0.0302 (0.0286-0.0350) 


15.99 


5.09 


3.97 


9408.24 


2995.26 


2334.04 


2.82 


0.28 


0.29 
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Table 4. Same as Table[3]but for method #2. 



Mag range 


Nb.obj 


Mass range 


Mid-mass 


dN 


errH 


errL 


dN/dM 


errH 


errL 


dN/dlogM 


errH 


errL 


12.0-12.5 


11 


u.ozuu- 


-U. J4UU 


0.5800 (0.5550-0.5890) 


11.00 


4.43 


3.28 


137.50 


55.35 


40.98 


2.26 


0.34 


0.35 


12.5-13.0 


54 


O ^/1AO 
U.J4UU- 


-u.4oyu 


0.5045 (0.4850-0.5125) 


54.00 


8.40 


7.33 


760.56 


118.30 


103.26 


2.95 


0.14 


0.15 


13.0-13.5 


91 


n A Aon 


n icon 
-u. joyu 


0.4290 (0.4025-0.4380) 


91.00 


10.58 


9.53 


1137.50 


132.23 


119.08 


3.05 


0.11 


0.11 


13.5-13.0 


106 


n is on 


-U. JUJU 


0.3460(0.3145-0.3595) 


106.00 


11.33 


10.28 


1232.56 


131.77 


119.58 


2.99 


0.10 


0.10 


14.0-14.5 


162 


U.jUjU- 


c\ linn 
-U.23UU 


0.2665 (0.2365-0.2800) 


162.00 


13.76 


12.72 


2219.18 


188.46 


174.22 


3.13 


0.08 


0.08 


14.5-15.0 


187 


0.2300- 


-0.1740 


0.2020 (0.1775-0.2140) 


187.00 


14.70 


13.67 


3339.29 


262.54 


244.03 


3.19 


0.08 


0.08 


15.0-15.5 


160 


0.1740- 


-0.1350 


0.1545 (0.1365-0.1655) 


160.00 


13.68 


12.64 


4102.56 


350.74 


324.08 


3.16 


0.08 


0.08 


15.5-16.0 


97 


0.1350- 


-0.1070 


0.1210 (0.1073-0.1305) 


97.00 


10.89 


9.84 


3464.29 


388.82 


351.29 


2.98 


0.1 1 


0.11 


16.0-16.5 


69 


0.1070- 


-0.0862 


0.0966(0.0858-0.1048) 


69.00 


9.35 


8.29 


3317.31 


449.60 


398.63 


2.87 


0.13 


0.13 


16.5-17.0 


58 


0.0862- 


-0.0710 


0.0786 (0.0671-0.0855) 


58.00 


8.66 


7.60 


3815.79 


570.06 


499.96 


2.84 


0.14 


0.14 


17.0-17.5 


46 


0.0710- 


-0.0549 


0.0630 (0.0534-0.0711) 


46.00 


7.84 


6.76 


2857.14 


486.79 


420.12 


2.61 


0.16 


0.16 


17.5-18.0 


20 


0.0549- 


-0.0486 


0.0517 (0.0463-0.0589) 


20.00 


5.56 


4.44 


3174.60 


881.78 


705.41 


2.58 


0.25 


0.25 


18.0-18.5 


18 


0.0486- 


-0.0440 


0.0463 (0.0414-0.0506) 


18.00 


5.33 


4.21 


3913.04 


1158.72 


915.89 


2.62 


0.26 


0.27 


18.5-19.0 


18 


0.0440- 


-0.0396 


0.0418 (0.0383-0.0461) 


18.00 


5.33 


4.21 


4090.91 


1211.39 


957.52 


2.59 


0.26 


0.27 


19.0-19.5 


9 


0.0396- 


-0.0368 


0.0382 (0.0356-0.0420) 


9.00 


4.12 


2.96 


3214.29 


1472.32 


1056.44 


2.45 


0.38 


0.40 


19.5-20.0 


15 


0.0368- 


-0.0339 


0.0353 (0.0330-0.0390) 


15.00 


4.97 


3.84 


5172.41 


1713.32 


1324.34 


2.62 


0.29 


0.30 


20.0-20.5 


22 


0.0339- 


-0.0311 


0.0325 (0.0305-0.0370) 


22.00 


5.77 


4.66 


7857.14 


2060.61 


1665.60 


2.77 


0.23 


0.24 


20.5-21.0 


4 


0.0311- 


-0.0294 


0.0302 (0.0286-0.0350) 


4.00 


3.18 


1.94 


2352.94 


1870.26 


1139.11 


2.21 


0.58 


0.66 



0.4 - 



0.2 - 



J: 

£ 0.0 



-0.2 



-0.4 




11 
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K1 
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WFCAM 



16 
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Figure 8. Difference in the K magnitude (K 1-K2) as a function of the K 1 
magnitude for all Pleiades member candidates with membership probability 
greater than 60%. The YJHK and JHK-only detections have been added 
too. Typical error bars on the colour shown as vertical dotted lines are added 
at the top of the plot. 



sample extracted by method #2 contains 1147 Pleiades member 
candidates, divided up into 978 stars (83.3%) and 169 brown 
dwarfs (14.7%). Similar percentages are obtained considering the 
sample of 947 high probability members (p^60%) identified in 
the probabilistic approach with 10.3% of brown dwarfs. Hence, 
the star (-0.6-0.08 M ) to brown dwarf (0.08-0.03 M ) ratio in 
the Pleiades lies between 5.4-6.3 and 8.6-9.3 if we consider a 3cr 
limit in the distance of the Pleiades (1 14.5-125.9 pc; lvan Leeuwenl 
2009). These numbers are in agreement with measurements (with 
slightly different mass ranges depen ding on the survey) derived 
from the field mass function ( 1 .7-5.3; lKroupal2002l;lchabrierl2005l ; 



lAndersen et al.ll20"06l) . young star-forming regions (3.0-6.4 for the 
Trapezium Cluster; 3.8-4.3 for a Orionis; 3.8 for Chamaeleon; 
8.3-11.6 for IC 348; lHiilenbrand & Can>entenl200(jt iMuench et"ai1 
20021; iLuhman et al.ll2003l ; lAndersen et alj|2006t iLuhmarj [2007k 



Lodieu et al. 2009), open c lusters (3.7 for the Pleiades and 4 .5 for 



M35; iBouvier et al. 1 1 19981; I Barrado y Navascues et al. 



2001), and 



Bate 2009, 



hydrodynamical simulations of star clusters (3.8-5.0; 
l201ll) . 



7.2 The cluster mass function 

In this section we adopt the logarith mic form of the Initial Mass 
Function as originally proposed by Salpeterl dl955l) : £(l°gio m ) 
= dn/dlog 10 (m) oc m~ a . We converted the luminosity in to a 
mass function using the NextGen models dBaraffe et al.ll 19981) for 
stars and brown dwa rfs more massive th an 50 Mj up (T e fi) and 
the DUSTY models IChabrier et alB oOO) for less massive brown 
dwarfs. The Z = 12-21.5 mag range translates into masses be- 
tween 0.62 and 0.03 Mq, assuming a revised distance of 120.2 pc 
( Ivan Le euwen 2009) and an age of 120 Myr for which the models 
are computed. 

We included in Fig.|9]errors in both the x-axis (logM) and y- 
axis (dN/dlogM) as follows. For the error bars on the masses, we 
considered three tim es the uncertainties on the age (125±8 Myr; 
Stauff er et al. and distance (120.2±1.9 pc; Ivan Leeuwenl 



20091) of the Pleiades given us a validity range of 3<r on the x-axis. 



Hence, we computed the masses with the 100 Myr NextGen and 
DUSTY isochrones shifted at a distance of 114.5 pc to define the 
lower limit and repeated the procedure with the 150 Myr isochrones 
for a distance of 125.9 pc as upper limits. (The uncertainties on 
the y-axis i.e. the dN/dlogM values are simply Gehrels error bars). 
The highest-mass point is incomplete due to the saturation of the 
GCS as are the two lowest-mass points of the mass function be- 
cause the last two magnitude bins in the Z-band are more con- 
taminated. Using these upper and lower bounds for the predicted 
masses for cluster members we refit the log-normal MF to exam- 
ine the effects o n the parameters of the fit (Fig.[l0]) Following the 
IChabrierj I I2003B Gaussian parameterisation (his Equation 12) 4e 
find a characteristic mass Tfic — 0.24^q Mq and a mass disper- 
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Figure 10. Mass function comparison between the Pleiades (shaded re- 
gion showing a lcr confidence interval from a weighted least-squares log- 
normal fit to the data presented here and extrapolat ed to lowe r masses) along 
with data from a Orionis (filled circles: [ L odieu et al.l2009l) and Upper Sco 
(crosses joined by a dashed line; lLodieu et alj|201ll) . These have been nor- 
malised at the peak of the MF. 



c-plot i 

feoosh ) 



sion a = 0.44 ± 0.01 dChabrierl d2005l) quotes m c = .25 and 
a = 0.55 for the disk system MF whereas Ichabrierl d2003l) quotes 
m c = 0.22 and a = 0.57). 

Overall we find that our Pleiades mass function is well rep- 
resented by a log-normal form over the 0.6-0.03 Mq mass range 
with a characteristic mass of 0.24 Mg (Fig.llOt. T his result is in 
agreement with all previous studies in the Pleiades l Moraux et al.l 



l200ltlTei et alj2002l ; lDobbie et al.l2002l:lDeacon & Hamblvl2004 ) 

over the same mass range and c onsistent with t he extrapolation 
of the system field mass function dChabrieri r2005) which can also 
reproduce preliminary densities of field L and T dwarfs found in 
large-scale surve' / s dMetchev et al.ll2008l: IBurningham e7ai1l20ld : 



lRevleetal.ll20ld ; iKirkpatrick et alj|201lh as displayed in Fig. [9] 



All these determinations of the mass functio n support the univer- 



excep t for the case of Taurus dBricefio et al. 200 2; Luhman et al. 

The latest hydrodynamical simulation of Bate (2011) is 
able to reproduce the observed field mass function (Kroupa 2002; 



excep t 
2003). 



sality of the IMF as discussed i n the review o flBastian et alJd2010h 

; latest hydrodynamical simulati' 
reprodu ce the observed field mass fu 
Ichabrierl 120051) with high confidence after inclusion of radiative 
feedback, in agreement with independent calcu lations dOffner et al.l 
120091 : lUrban et aljioiol : iKrumholz etafll201ll) . 

It is interesting to compare an extrapolation to lower masses 
of the log-normal Pleiades mass function derived here with results 
from other younger GCS targets (Fig. |10| . While the Pleiades is 
measured in a higher mass range and shows a closely log-normal 
form, the I MF from a carefully cleaned spectroscopic sample in 
Upper Sco dLodieu et al.l 1201 lh penetrates to lower masses and 
seems to be much shallower in the substellar regime, implying 
the presence of more brown dw arfs. Results for the sigma Orio- 
nis cluster show a similar trend dLodieu et al.ll2009l) . although to 
a less obvious extent. Of course, these mass functions have been 
derived assuming that the evolutionary models accurately predict 
colours, bolometric magnitudes and temperatures at the different 
ages (especially ages < 10 Myr) and assuming also that any sys- 
tematic errors introduced by not accounting for unresolved binarity 
cancel in the comparison. In Fig. [TT]we show an 'a;-plot' equiva- 
lent to the 'alpha plot' with x = a — 1), i.e. a plot of the gradi- 
ent of mass function as a function of mass, in order to compare 



equivalent to the 'Alpha plot' (where x = a — 1, 
2005)) showing the gradient of the mass function as a 
function of mass. The straight line is the gradient of the log-normal fit 
to the Pleiades data presented herein, while th e filled circle is th e re- 
cent GCS results from a Orionis (a=0.5±0.2; lLodie u et al. 200j^), the 
cross is from a po wer-law fit to the most recent Upper Sco results 
(x = -0.6 ± 0.08; lLodieu et alj|201ll), the open dia mond for Alpha Per 
n Ban 



(a=0.59±0.05) from Ba rrado v Navascues etall (l2002), and the filled trian- 



iKraus & Hille nbrand 2o67r. lBoudreault et alJl2010tlBaker etai]|2O10l) . The 
arrow indicates the upper limit to the mass function of the field brown dwarf 
population IBurningham et alj2010l) . Horizontal dotted lines indicate mass 
ranges for power-law IMF determinations. 



with the constraints set on the field IMF gradient by the UKIDSS 
Large Area Survey T-type b rown dwarf searches (e.g . lLodieu et al.l 
l2007l : IPinfield et all l2008t IBurningham et alJ l20ld) ). While the 
mass functions of the more aged populations like Praesepe 
dKraus & Hillenbrand 120071: iBoiidreault et ail l20ld: iBaker et all 



20ld) or the field dMetchev et alj|2008t IBurningham et al.l l2010l: 
Revle et al.|[2oTcl ; IKirkp atrick et al.ll201ll) show gradients consis- 



tent with the Pleiades log-normal, power-law fits to the very young 
Upper Sco (x = -0.4 ± 0.08;lLodieu et al.ll20TTh and a Orionis 
(x=-0.5±0.2; lLodieu et al.ll2009h clusters" are flatter. This may be 
evidence of a variation in the substellar IMF, or simply an artefact 
of systematic errors in the evolutionary model predictions at very 
young ages or age spreads. We should emphasise that the upper 
limit set for the field mass function by IBurningham et alJ d2010h 
remains under debate (due to the lack of constraints on ma sses 
and ages) as the CFHT brown dwarf survey dRevle et al .120 id) and 
the p reliminary densities determined by WISE (K irkpatrick et al.l 
l201ll) suggest p ositive a values for th e mass function in the T dwarf 
regime (see also lMetchev et al]|2008l) . 



8 SUMMARY 

We have presented the outcome of a wide (~80 square degrees) and 
deep (J ~ 18.8 mag) survey in the Pleiades open cluster as part of 
the UKIDSS Galactic Clusters Survey Data Release 9. The main 
results of our analysis can be summarised as follows: 

• we recovered Pleiades member candidates previously pub- 
lished and updated their membership assignations 

• we selected photometrically and astrometrically potential 
Pleiades member candidates using two independent but comple- 
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mentary methods: the probabilistic analysis and a more standard 
method combining photometry and proper motion cuts 

• we derived a BD binary fraction around 25.6±4.5% in the 
0.075-0.03 Mq mass range with a difference of a factor of two be- 
tween high-mass (0.075-0.05 M ) and low-mass (0.05-0.03 M Q ) 
brown dwarfs for separations less than ~ 100-200 au 

• we investigated the if-band variability of Pleiades members 
and found virtually no variability at the level of 0.08 mag 

• we derived the luminosity function from both selection meth- 
ods and found no difference within the error bars 

• we derived the Pleiades mass function which is best fit by a 
log-normal function peaking at 0.16-0.20 M Q in the 0.6-0.03 M 
range, in agreement with previous studies in the cluster and the 
extrapolation of the system field mass function 

This paper represents a significant improvement in our knowl- 
edge of the Pleiades population and the cluster mass function in 
the substellar regime over the full cluster. We believe that this pa- 
per will represent a reference for many more years to come. We 
will now extend this study to other regions surveyed by the GCS 
to adress the question of the universality of the mass function us- 
ing an homogeneous set of photometric and astrometric data. Fu- 
ture work to constraint current models of star formation include the 
search for companions to investigate their multiplicity properties, 
the determination of the radial velocities of Pleiades members, and 
deeper surveys to test the theory of the fragmentation limit. 
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Table Al. Sample of 1379 known Pleiades member candidates previously published in the literature and recovered in GCS DR9. We list the equatorial 
coordinates (J2000), GCS ZY J H K1K2 photometry, proper motions (mas/yr) and their errors, reduced chi-squared statistic of the astrometric fit for each 
source (x 2 value), membership probabilities, and names from the literature. Note that the 312 sources without membership probabilities are divided into four 
groups: 190 non members (NM) detected in ZYJHK but not satisfying our photometric and astrometric criteria, 42 objects without Z + Y photometry 
(noZY), 74 objects without Z only (noZ), and 6 sources without Y only (noY). Pleiades member candidates are ordered by increasing Right Ascension. This 
table is available electronically in the online version of the journal. 



R.A. 


Dec. 


Z 


Y 


J 


H 


K\ 


K2 


fi a cosS 1 err 


/is ±err 


x 2 


Prob 


Name 


03 27 54.26 
03 29 58.76 


+24 56 10.9 
+23 22 18.3 


13.808 
13.640 


13.400 
13.198 


12.820 
12.672 


12.255 
12.244 


12.005 
11.843 


99.999 
11.851 


16.62±6.95 
21.18+3.41 


-43.60+6.95 
-38.83+3.41 


0.47 
0.47 


0.93 
0.81 


DH003 
DH009 


04 05 13.75 
04 06 29.99 


+24 08 42.7 
+22 33 43.6 


14.983 
14.376 


14.471 
13.856 


13.846 
13.201 


13.261 
12.623 


12.933 
12.273 


12.917 
12.269 


19.77±3.38 
14.83±5.07 


-41.50+3.38 
-33.96±5.07 


0.54 
0.55 


0.94 
0.14 


DH915 
DH2004.916 



Table Bl. Coordinates (J2000) and names of 544 previously-known Pleiades member candidates not recovered in the GCS DR9. Pleiades candidates are 
ordered by increasing right ascension. This table is available electronically in the online version of the journal. 



R.A. 


Dec. 


Old Names 


03:27:42.06 


+23:48:13.3 


PELS 121 


03:28:01.56 


+23:04:42.6 


DH004 


04:05:09.44 


+23:28:59.0 


DH913 


04:05:13.72 


+22:18:19.0 


DH914 



Table CI. Coordinates (J2000), near-infrared (ZY JHKXKT) photometry with the error bars, proper motions with errors for all new Pleiades member can- 
didates identified in the UKIDSS GCS DR9 with the probabilistic and standard selection methods. The penultimate column gives the membership probability 
if the object was selected with the probabilistic method. The last column lists the source (S) of the object: "1" means that the object was identified in the 
probabilistic approach, "2" means that the object was identified with method #2, "12" means that the candidate is common to both selection method, and 
"3" means that this is a faint YJHK or JHK candidate. Pleiades member candidates are ordered by increasing Right Ascension. This table is available 
electronically in the online version of the journal. 



R.A. Dec. 


Z 


Y 


J 


H 


Kl 


K2 


// a cosS 




Prob 


S 


03 24 59.74 +25 34 04.5 


16.878+ 0.009 


16.243+ 0.007 


15.582+ 0.007 


14.987+ 0.010 


14.606+ 0.005 


99.999+99.999 


38.59+7.21 


-49.57+7.21 




2 


03 25 30.92 +24 51 39.9 


20.316+ 0.126 


19.478+ 0.096 


18.600+ 0.079 


18.210+ 0.160 


18.170+ 0.120 


99.999+99.999 


10.72+31.33 


-45.76+31.33 


0.60 


1 


04 10 54.54 +26 01 42.4 


19.642+ 0.091 


99.999+99.999 


17.748+ 0.036 


17.002+ 0.027 


16.417+ 0.022 


99.999+99.999 


17.38+5.99 


-47.11+5.99 


0.66 


1 


04 11 03.84 +23 15 48.9 


17.766+ 0.017 


17.254+ 0.014 


16.599+ 0.014 


15.986+ 0.017 


15.599+ 0.021 


15.623+ 0.012 


18.36+6.81 


-46.34+6.81 


0.61 


1 
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Table Dl. Coordinates (J2000), near-infrared (ZYJHK1K2) photometry, and proper motions (in mas/yr) for substellar multiple system candidates identified 
photometrically in the Pleiades cluster 



R.A. 



Dec. 



Y 



H 



K\ 



K2 



H a cos8 



03 25 38.73 
03 31 20.71 
03 32 11.55 
03 33 49.22 
03 34 38.61 
03 36 01.95 
03 36 03.85 
03 36 53.30 
03 40 11.98 
03 40 45.16 
03 40 53.66 
03 41 40.91 
03 41 42.41 
03 41 54.16 
03 43 34.49 
03 44 14.65 
03 44 35.16 
03 44 35.90 
03 45 09.46 
03 45 31.37 
03 45 37.76 
03 45 41.27 
03 45 50.66 
03 46 05.11 
03 46 15.11 
03 46 20.27 
03 46 22.25 
03 46 26.09 
03 46 27.10 
03 47 11.79 
03 47 20.48 
03 48 04.67 
03 48 31.53 
03 48 35.20 
03 48 50.45 
03 48 57.41 
03 50 52.17 
03 51 38.96 
03 53 55.13 
03 55 27.06 
03 56 52.31 
03 58 00.62 
03 58 17.43 

03 59 59.85 

04 00 03.21 
04 00 08.16 
04 00 50.52 
04 01 28.43 
04 01 39.83 
04 01 50.95 
04 09 17.80 



+22 57 39.8 
+25 57 33.6 
+21 27 55.7 
+19 59 52.0 
+24 51 28.5 
+27 1 1 04.7 
+22 52 03.5 
+26 34 27.9 
+21 48 31.8 
+27 50 40.4 
+28 21 11.5 
+25 54 24.1 
+23 54 57.1 
+23 05 04.7 
+25 57 30.6 
+23 49 40.0 
+25 13 42.8 
+23 34 41.9 
+23 58 44.7 
+24 52 47.4 
+23 43 50.1 
+23 54 09.7 
+24 09 03.5 
+23 45 34.9 
+26 46 48.8 
+23 58 18.9 
+23 52 26.6 
+24 05 09.5 
+21 48 22.6 
+24 13 31.3 
+19 54 25.5 
+23 39 30.1 
+24 34 37.2 
+22 53 42.1 
+22 44 29.8 
+23 13 59.1 
+23 27 11.2 
+24 30 44.8 
+23 23 36.1 
+25 14 45.8 
+25 10 05.1 
+21 18 20.8 
+22 1 1 52.7 
+25 08 53.6 
+22 24 46.0 
+22 32 01.1 
+23 43 52.9 
+23 30 59.6 
+22 47 53.7 
+22 59 15.5 
+26 03 31.2 



15.890 
16.847 
16.385 
18.539 
17.097 
16.248 
18.072 
17.046 
16.578 
17.074 
18.880 
16.893 
16.171 
17.349 
16.571 
15.892 
17.656 
16.307 
16.974 
17.332 
16.240 
17.166 
17.478 
17.229 
16.624 
19.259 
17.120 
16.810 
19.797 
16.305 
17.011 
17.014 
19.218 
16.176 
16.562 
16.835 
17.690 
18.711 
16.947 
16.118 
16.146 
17.399 
16.259 
16.368 
16.534 
16.449 
16.184 
16.318 
16.669 
16.363 
17.120 



15.269 
16.068 
15.652 
18.104 
16.247 
15.492 
16.873 
16.171 
15.884 
16.218 
17.570 
16.001 
15.464 
16.376 
15.727 
15.245 
16.584 
15.672 
16.250 
16.330 
15.456 
16.189 
16.582 
16.347 
15.838 
18.174 
16.323 
15.967 
18.650 
15.554 
16.054 
16.054 
17.883 
15.452 
15.825 
16.033 
16.638 
17.407 
16.027 
15.402 
15.491 
16.380 
15.503 
15.695 
15.875 
15.820 
15.380 
15.539 
15.899 
15.631 
16.480 



14.601 
15.309 
14.906 
16.453 
15.459 
14.774 
15.913 
15.379 
15.169 
15.404 
16.490 
15.180 
14.709 
15.522 
14.909 
14.547 
15.662 
14.985 
15.438 
15.465 
14.715 
15.360 
15.705 
15.522 
15.032 
17.034 
15.518 
15.160 
17.374 
14.792 
15.210 
15.283 
16.715 
14.745 
15.098 
15.222 
15.700 
16.400 
15.172 
14.649 
14.756 
15.545 
14.768 
15.004 
15.128 
15.102 
14.647 
14.772 
15.153 
14.906 
15.759 



13.974 
14.745 
14.339 
15.767 
14.920 
14.234 
15.240 
14.804 
14.580 
14.816 
15.852 
14.574 
14.124 
14.975 
14.359 
13.957 
14.985 
14.376 
14.872 
14.839 
14.172 
14.782 
15.095 
14.851 
14.495 
16.269 
14.917 
14.584 
16.564 
14.261 
14.615 
14.704 
15.977 
14.185 
14.531 
14.673 
15.060 
15.718 
14.569 
14.071 
14.192 
14.931 
14.193 
14.433 
14.550 
14.503 
14.089 
14.202 
14.575 
14.311 
15.003 



13.618 
14.321 
13.949 
15.433 
14.445 
13.810 
14.679 
14.355 
14.187 
14.321 
15.194 
14.122 
13.686 
14.415 
13.909 
13.546 
14.448 
13.990 
14.424 
14.354 
13.756 
14.305 
14.580 
14.385 
14.037 
15.650 
14.474 
14.118 
15.848 
13.841 
14.152 
14.256 
15.357 
13.770 
14.141 
14.194 
14.505 
15.168 
14.088 
13.671 
13.771 
14.445 
13.787 
14.029 
14.132 
14.098 
13.660 
13.788 
14.169 
13.912 
14.603 



13.628 
14.328 
13.920 
15.492 
14.466 
13.782 
14.679 
14.353 
14.178 
14.333 
15.199 
14.125 
99.999 
14.418 
13.901 
99.999 
14.460 
13.985 
14.410 
14.326 
13.742 
14.309 
14.560 
14.404 
14.088 
15.585 
14.472 
14.098 
15.925 
13.850 
14.140 
14.238 
15.312 
13.772 
14.143 
14.194 
14.555 
15.122 
14.081 
13.650 
13.801 
14.423 
13.778 
14.075 
14.141 
14.104 
99.999 
13.788 
14.157 
13.906 
99.999 



24.45±5.11 
21.73±3.39 
13.64±3.89 

I. 67±6.44 
13.39±2.69 
15.79±3.79 
22.46±3.15 
15.54±3.33 
24.16±2.54 
12.37±3.34 
13.91±4.06 
16.94±2.26 
14.07±2.31 
18.19±2.30 
20.92±2.25 
14.46±2.31 
19.33±2.26 
16.94±2.23 
16.07±2.25 
16.69±2.24 
20.91±2.23 
17.46±2.24 
16.01±2.26 
15.88±2.25 
2().73±2.97 
12.97±2.40 
17.65±2.25 
19.41±2.24 
20.95±2.98 
16.88±2.23 
28.11±5.10 
16.07±2.24 

II. 92±2.37 
15.22±2.15 
15.64±2.16 
14.93±2.17 
19.89±2.21 
23.01±2.34 
19.17±2.28 
15.24±2.24 
16.66±2.50 
24.81±3.43 
20.64±2.89 
14.87±2.51 
15.92±2.87 
16.31±2.88 
24.52±3.56 
24.55±3.37 
21.16±3.39 
20.72±3.38 
6.56±4.49 



-47.02±5.11 
-38.42±3.39 
-40.00±3.89 
-40.37±6.44 
-37.21±2.69 
-42.70±3.79 
-46.19±3.15 
-34.04±3.33 
-40.67±2.54 
-40.37±3.34 
-31.57±4.06 
-42.13±2.26 
-47.37±2.31 
-44.74±2.30 
-47.72±2.25 
-40.15±2.31 
-44.97±2.26 
-44.39±2.23 
-42.23±2.25 
-40.30±2.24 
-45.45±2.23 
-44.47±2.24 
-40.58±2.26 
-39.21±2.25 
-42.79±2.97 
-35.00±2.40 
-38.04±2.25 
-38.71±2.24 
-48.67±2.98 
-41.27±2.23 
-39.39±5.10 
-44.27±2.24 
-46.68±2.37 
-45.62±2.15 
-42.06±2.16 
-36.62±2.17 
-41.33±2.21 
-40.47±2.34 
-44.72±2.28 
-39.66±2.24 
-38.22±2.50 
-32.13±3.43 
-34.95±2.89 
-35.69±2.51 
-33.98±2.87 
-39.63±2.88 
-40.90±3.56 
-39.20±3.37 
-42.62±3.39 
-38.98±3.38 
-33.73±4.49 
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